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Abstract 

We extend the results of a previous paper on the Gross-Pitaevskii description of rotating Bose- 
Einstein condensates in two-dimensional traps to confining potentials of the form V{r) = r" , 2 < s < 
oo. Writing the coupling constant as we study the limit e ^ 0. We derive rigorously the leading 
, asymptotics of the ground state energy and the density profile when the rotation velocity fl tends to 

' infinity as a power of 1/e. The case of asymptotically homogeneous potentials is also discussed. 
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> ■ 1 Introduction 

00 

In a previous investigation of rapidly rotating Bose-Einstein condensates in two-dimensional anharmonic 
traps [CDY| we considered the case of a 'flat' trap with a rigid boundary confining the condensate to a 
disk of finite radius. The present paper is a sequel to this work, extending the results to homogeneous 
I trap potentials of the form V(r) = with 2 < s < oo. The flat trap corresponds to the limiting case 

s — > oo. For s < oo the system is no longer confined to a bounded region when the coupling constant 
I tends to infinity and a suitable scaling of the variables is necessary to obtain a well-defined limit. This 

gives rise to additional features and requires some modifications that are dealt with in the present paper. 
As in [CDY] we identify three different parameter regimes depending on the way the rotation velocity is 
^ ' scaled with the interaction. 

^ I The present study is carried out strictly within the Gross-Pitaevskii (GP) framework in contrast to 

the recent paper [BCPYj where the main emphasis is on many-body aspects and the GP description is an 
auxiliary tool. The trap potentials considered there are three dimensional and not necessarily rotationally 
symmetric. A two-dimensional potential of the form is an instructive special case that can be analyzed 
in more detail than the general case. 

As in [CDY], where a discussion of the general context and an extensive list of references can be 
found, the starting point is the two-dimensional Gross-Pitaevskii energy functional for the wave function 
^! of the condensate which in our case can be written as 

£GP[^]= f df(|V*|2 + r^|5'|2-0(£)^*L* + ^l. (1.1) 
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Here L is the third component of the angular momentum (i.e. L = —id/d-d in polar coordinates (r, i?)), 
the angular velocity and e is a non- negative, small parameter. The wave function is normalized so 
that l^'P ~ 1- Units have been chosen so that h = 2m = 1, where m is the particle mass, and such 
that the coefficient in front of r'^ is simply 1 . 

If s < oo the condensate spreads out indefinitely in the Thomas-Fermi (TF) limit e and the 
density j^'p tends uniformly to zero. Non-trivial results can be obtained, however, by rescaling all 
lengths by an e-dependent factor. We write r = kr' , ^'(r) = '$'{r')/k with k = e~2/(s+2) ^^^^ define 

£^'^[-i'']=e~^£^'^[^]. (1.2) 
The functional £'~'^ is the one we shall study. Dropping the primes of the arguments it is explicitly given 

by 

£GP[^] ^ J ^ ||v^|2 _ + + |^|2)| ^ (13) 

where the scaled angular velocity is given by 

w(e) = e-^O(e). (1.4) 
The functional (|1.3p is defined on the domain 

pGP = (= H^R^) I r^|*|2 G l1(M2)}. 

We set 

£;?P= min S'^^m (1.5) 

*eP<5P,||*||2=i 

and denote by ^Pp^ a corresponding minimizer, which may not be unique [S]. 

In the following we study the leading order asymptotics of the ground state energy and density for 
the functional £^^[5'] as e ^ 0. The ground state behavior for the original GP functional (|l.ip follows 
by scaling: If we set E^^ = inf £'^^[5'] and denote by '^f^ any ground state of (jl.ip . one has 
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e^Ef^ *?^(^) = e^^f^ (^e^r^ . (1.6) 



Note that the 'flat' trap case studied in FB] and [CDY| can be formally obtained by taking the limit 
s ^ oo of (|1.3p : In this limit the scaling factor e~2/(s-i-2) converges to 1, the external potential to oo 
for r > 1 and to for r < 1, and the rescaled angular velocity uj{e) to ri(e). The formal limit s — > oo 
corresponds to a 'flat' trap with Dirichlet conditions at the boundary rather than the Neumann conditions 
considered in [CDY| . As noted in [CDY| both boundary conditions give the same results in the TF limit 
e^O. 

As in [CDY| we rewrite the GP functional in the form 

8--m - I df {I (V - ^A^) ^\\^^ir^ + m - ^^ifl!:!^ | (1.7) 
with the vector potential 

A,ir) EE ^^e- X r (1.8) 

where is the unit vector in z-direction. The behavior of the GP functional as e — s- depends on the 
way the angular velocity u) scales as a function of e. We distinguish three cases: uj <^ 1/e, oj ^ 1/e 
and Lu ^ 1/e. It is convenient to write lu = loq/s] the three cases then correspond to ^ 1, ^^o ~ 1 
and luq ^ 1. In the next section we discuss for fixed ujq the TF functional that is obtained from the 
GP functional by dropping the flrst (kinetic) term in (|1.7D . This functional and its limits for cjq ~> 
and LUo —^ 00 describe the asymptotics of (jl.7p for e — > as summarized in Section 3. In Section 4.1 we 
present the proofs for the regimes <C 1 and uj ^ 1/e and in Section 4.2 for oj ^ 1/e. 

Since the length scale e~^/^''+^^ tends to infinity in the TF limit, it is clear that for leading order 
calculations only the asymptotic behavior of the confining trap potential for large arguments matters. 



In the last section we indicate how our proofs can be extended to include potentials V{r) that are 
asymptotically homogeneous in the sense of [LSYlj . i.e., such that for A > 1 

\X-''V{Xr)-r'\<cX-''{l + r'') (1.9) 

for some constants k, c > 0, uniformly in r G M^. 

2 The TF functional and its properties 

The TF functional depends on the density alone and is for fixed loq defined as 



on the domain 

= {pe L\R^) \p>0,r'pe L1(R2)} . (2.2) 
By standard methods there is a unique minimizer 



p-(.)4 



(2.3) 
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where the chemical potential pT'^ is fixed by = 1. The ground state energy associated with (|2.3p 

is 

i?TF= inf £:TF[ ] ^^TFj TF]^ (2.4) 

pei3^^.||p||i=i 

Since — r*+Wgr^/4 — oo as r — » oo, due to the condition s > 2, the minimizer (j2.3p is always compactly 
supported, i.e., supp(p'^^) C Sr, for some R < oo depending on ujq where Br denotes a two-dimensional 
ball centered at the origin with radius R. 

As in the case of a flat trap the density p"""^ develops a 'hole' (a disk centered at the origin where the 
density (j2.3p vanishes) when luq exceeds a certain critical value ujq.c- Both the outer radius i?out of the 
support of p"^^ and the inner radius i?in (if present) increase with luq as discussed below. In the flat trap 
only i?in increases with luq while -Rout is fixed from the outset. 

The chemical potential depends also on uiq and because — r" + cjgr^/4 is monotone increasing in 
ujQ, it is clear that p"^^ is monotone decreasing in uiq (see also (|2.12p and (|2.14p ). For small ujq we have 
/i'^F > 0, whereas vanishes and changes sign at ujq = ^o,c- 



2.1 Support of p for ujq = const. 

In order to study the support of p"^^ it is convenient to consider the function f{z) = p^^ — z*/^ + (cjq/4)z 
where z = > 0, so that p^^ = f{r^)+. Since d^//dz^ — — s(s — 2)zt~/4 < for z e (0,cx)), the 
function / is strictly concave for z > 0. Moreover /(O) = p'^^ and lim^^oo f{z) = — oo. Hence, if 
p'^^ > 0, i.e., utQ < ujQ.c, then /(O) > and there exists a unique Zout > such that /(zout) = 0. The 
support of p^^ in the radial coordinate is the interval [O,i?out], with i?out = V^out- In the opposite case 
p^^ < 0, i.e., Wo > ujo.c, we have /(O) < but sup / > (since / p'^^ > 0) and concavity implies the 
existence of two positive solutions of /(z) = 0. The support of p"^^ in the radial coordinate is then an 
interval [i?in,-Rout] for some < Rin < Rout- Note also that 

0</'(i?L) = -^i?f-2 + ^| and 0>/'(i?L) = -^i?out + ^- (2.5) 
From the L"'^— normalization of p"^^ and 



, ,2 r2 
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(2.6) 



we get for ujo < wq,, 



T>s+2 , ,2 p4 -1 

2{s + 2) 16 ~ tt' ^ ■ ^ 

whereas for ojq > tt>o,c we have 

, ,2 p2 , ,2 n2 

A* =-Kin = ^out ; l^-fSj 

2(s + 2) 16 tt' ^ ' 

The radii i?out and i?in are determined by solving (|2.7p . or (|2.8p together with (|2.9p . While explicit 
formulas can in general not be given, the outer radius i?out.c for m"'"^ — and the critical angular velocity 
wOjC for the creation of a hole are easily obtained from ()2.6p and ()2.7p : 

i?out,c=(^^)'^ (2.10) 

with 



Wo,, 



4(s + 2) 



(2.11) 



_7r(s-2)_ 

In the 'flat' trap case, i.e., for s oo, the expression for cjo.c simplifies to 4/^/?? as in ICDYj . 

We now show that fi^^ is monotonically decreasing and i?in and i?out monotonically increasing as ljq 
increases. Consider first the case ujq < ojq^c, i.e., /i"^^ > 0. Differentiating the normalization equation 
J p^^ = 1 with respect to i = cjq/4, using (|2.6p . gives 

a/i^F/^i = < 0. (2.12) 

Differentiating (j2.7p gives 

9i?out/9i = (si?:-t' - 2t)-ii?out/2 (2.13) 
and hence dRout/dt > because of (|2.5p . 



In the case cjq > ^^o,c, we again differentiate the normalization condition for p^^ , this time using 
and obtain 

dp^^/dt = -{Rl^^ + i?fJ/2 < 0. (2.14) 
Moreover, by taking the derivative of p.Sp w.r.t. t, we have 

— V*^i„ ^I'J^in) ^in — l,*^out ^^^out; ^outJ 

which, combined with (|2.14p . yields the inequalities 

n p n p 

(si?f-i-2ti?i„)i^ <0 and (sK^^i - 2ii?out)^^ > 0. 

By (|2.5p this implies 

^>0 and ^>0. (2.15) 

Altogether, we have thus seen that p^^ has compact support contained in Bn^^^. If ujq < luq^c, the 
support coincides with Br^^^, while, for luq > wo,c, it is the annulus {r : Rm <r< i?out}. Both R^ and 
-Rout grow as luq increases. 



2.2 The quartic trap 

Let us consider for illustration the quartic trap V{r) = . The critical angular velocity and radius are 
given by 



Wo,, 



12 



1/6 



i2y/« 



The equations for the inner and outer radii i?in and i?out reduce to equations of third order for the squares 
of the radii. If ujq < cjq.c the chemical potential and the outer radius are 



TF _ J_ I i li^a) 
^ " 64 1 2 



1/3 



1/3 



64 



and 



Rnut — 



1/3 



1/3 



1/2 



with q{ujo) = 6144 + ttujq + 64^3^^3072 + ttwq. If ujo > a;o,c, we obtain 
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1 /12 



2/3 4 

64' 
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TT 



1/3 



1/3 



(Note that the two expressions for fi^^ and i?out are equal for ojq — ojq.c whereas -Rin = 0.) For luq > wg,, 
the TF minimizer is 



16 



r 



1 /12 



8 V TT 



2/3' 



2.3 Support of p^^ for cjq oo 

The radius at which the density is maximal can be explicitly calculated from 

2s 



(2.16) 



It is clear that < Ry^ < i?out and all radii tend to infinity if — > oo. We shall now show that 
^out — -Rin tends to zero in this limit. For s > 4 also R^^t — Rf^ tends to zero and the density therefore 
to infinity. 

It is convenient to scale all lengths by using as a unit, i.e, to write r = R^x. The scaled TF 
minimizer fp-^ (x) = Rf-^p^^ (r) is 



p^^ix) 



1 fc^l 



2 V 2s 



-TF s 
p — X 



sx 



(2.17) 



with the scaled chemical potential 



2s 



TF 



(2.18) 



We also denote x\n = Rin/ Rm and a;out = Rout /Run so that < Xi^i < 1 < a;out and the maximum of p^^ 
is attained at x = 1. 

In the same way as (|2.14p was derived we have 



dp 



r,TF 



2s 



(s + 2) 



7rs(s-2)(a;2^t-a;2j 



< 



(2.19) 



so that p. is a decreasing function of ujq. Moreover since p (1) > 0, one has the bound 



r.TF 



< - - 1. (2.20) 



Defining 

h{x) = x' -{s/2)x'^ + {s/2)-l, (2.21) 
the scaled radii Xm and Xout are determined by 

h{xin) = HXout) (2.22) 

together with the normahzation condition for fp-^ which expressed in terms of h is 

^{xout'~xl)~ £'^\ix)xdx = 7r-'(^^^ (2.23) 

The right hand side of (j2.23p tends to zero as ^ oo and the integral on the left hand side is always 
strictly less than the first term for Xin strictly less than Xout- Since Xin < 1 < Xout and h is continuous 
with h(x) = only for a; = 1 it is clear that Xm and Xout must both tend to 1 as o^o — > oo and by the 
normalization the density fp-^ approaches a delta function concentrated on the circle with radius 1 . Note 
also that 

/i(xi„) =/i(xout) = + ^ -1 (2.24) 

so that 

— (2.25) 

In order to estimate the rate of the convergence of the density to a delta function we make a Taylor 
expansion of h{x) around a; = 1: 

h{x) = is(s - 2)(1 - x)^ + 0(|1 - xf). (2.26) 

Writing Xin = 1 — (5 + o{6), Xout ^ 1 + S + o{S), where S is the deviation from 1 to leading order in the 
small parameter on the right hand side of ()2.23p . the normalization condition (I2.23P gives 

Multiplying 6 with ~ Wq^^*"^"* we see that 

(i?out - i?i„) ~ c^o^'^'^'^/'^'^'^ (2.28) 

Thus also the original density p^^ is supported on an annulus whose thickness tend to zero. The area of 
the support is 

7r(i?L - Rl) - (2.29) 
which increases for 2 < s < 4 but tends to zero for s > 4. 

2.4 TF energy asymptotics for a;o — >^ oo 

The scaled density p^^ is the minimizer of the scaled TF functional 



£'^^[p\ = I dx 



sx 



(2.30) 



with corresponding energy — £"^^[p^^]. As shown in the previous subsection p^^ converges to a 
delta function on the unit circle as t^o ~* oo. The behavior of the energy in this limit is given in the 
following proposition. 



so that 



Proposition 2.1 (TF energy for luq oo) 

For ujQ — !■ oo, 

i?TF = 1 - £ + O ^<=-^'j . (2.31) 

Proof: The lower bound is simply obtained by neglecting the positive last term in (|2.30p and using the 
inequality a;" — sx'^/2 > 1 — s/2. For the upper bound we use a trial function of the form 

p^{x)^r'jir\i-^')) (2.32) 

where j is a smooth non-negative function supported in [—1/2,1/2] satisfying the normalization 
TT J drj{r) — 1 and < ^ < 1. One can easily estimate ||p{||2 < CC~^, and exploiting the Taylor 
expansion of — sx- /2 around x — I we have 

dx(^x^-'-^y,<i-'i+c'e 

^~^^[/5d < 1 - i + + C" u^-'^r'- (2.33) 
Optimization with respect to the parameter ^ yields the desired upper bound. 

□ 

3 Main results 

3.1 The regime u <^ 1/e 

For CO <^ 1/e, the GP ground state energy and density are approximated to the leading order by the 
corresponding quantities in the non-rotating case, exactly as in Proposition 2.3 in [CDY) . The TF 
functional without rotation, i.e., for luq = 0, is given by 

Sj^ip]^ I dfp(r^+p). 

We denote by E'J^ its ground state energy and by 

pl^ir)^\[p^^^r^]^ 

the corresponding minimizer. 

Proposition 3.1 (GP energy and density asymptotics for ^ 1/e) 
For any w(e) such that lim£__,o £to{s) = and for e tending to zero, 

e^Er^Ej^ + o{l), (3.1) 

III*?'!' -P^'IL^CM^) -«(!)• (3-2) 

Proof: The lower bound for the ground state energy is actually trivial, since it is sufficient to neglect 
the first positive term in (|1.7p to obtain 

In order to get an appropriate upper bound we test the GP functional on the (real) GP minimizer for 
w = and obtain Ef^ < E^^\^^o. The result is then a consequence of Lemma 2.3 in [LSYlj and we get 
the bound (see Eq. (2.18) in [LSYlp 

The density convergence is a simple corollary (see the proof of Proposition 2.3 in |CDY| and Theorem 
2.1 in [LSTI]). 

□ 



3.2 The regime uj ^ 1/e 

We now assume that a;(e) = ^o/£ with > a finite constant. The analogs of Theorem 2.1 and 
Corollary 2.1 in [CDYj are the following: 

Theorem 3.1 (GP energy asymptotics for lo ^ 1/e) 
For any Wq > o.'^^d for e tending to zero, 



s'EY^ = E''' +0(e|loge|). 

Corollary 3.1 (GP density asymptotics for u! ~ 1/e) 
For any ojq > and for e tending to zero, 



(3.4) 



\\\^T?-p'''\\lh-^.)^0{^M^\). (3.5) 
The asymptotics of the energy and density for the original functional p.ip quantities is then given by 

2 



where the convergence of the density is in the norm topology of L^{M?). 

For s < oo, the condensate is not confined to a bounded region and is a function supported 

on the whole of M^. From Corollarv l3.1l it follows immediately that j^'^^p is small outside the support 
of p^^ in norm, i.e., 

df = 0(e|loge|) (3.6) 

2\supp(pTF) 

but much more can be shown, namely that |^^^| is pointwise exponentially small outside the support of 

TF 

Theorem 3.2 (Exponential smallness of the GP density, w ~ 1/e) 

For any > 0, r e 7^°^* = {r S | r > /?out + e^^^} and for e sufficiently small, 



j 



exp 



C' dist(f,9T™t)2 



r5/6 



(3.7) 



Furthermore, for any ujo > ujq^c the same estimate holds for f G T"^ = | r < i?in — e^^^} and 

Qj-out replaced with dV''. 



3.3 The regime cij ^ 1/e 

For convenience, in particular for the statement of Theorem 13.41 below, and comparison with |CDY] we 
assume that lo increases as a power of 1/e, i.e., that (jj(e) — wi/e^"''" with some constants U!i,a > 0. 
This means that we take utQ — eaj(e) — uji/e". Theorem l3. 31 holds true for general ujo{e) = ew(e) ^ oo if 
uji/e" is replaced by i^oi^)- 

In the regime w ^ 1/e the limiting functional is still given by (|2.ip . but since loq now depends on e 
this is also the case for the TF ground state energy and density. We thus use the notations Ej^ and 
p^^ . Proposition 12. II vields the ground state energy asymptotics for the functional Ej^ , i.e., 

.a.j..(|i)*(,_£),o(.«^)^ (") 

The following Theorem describes the GP ground state energy asymptotics. 

Theorem 3.3 (GP energy asymptotics for uj ^ 1/e) 
For any wi , a > and e tending to zero, 

eSl = (il) (l - i) + ^ (^'^) + ^ ^"^) • (3-9) 



The first one, of order e 3{s-2> ^ jg actually the expected 



Therefore, as long as a < ^^'^ ^-^ 



Note the occurrence of two remainders in 
optimal one, since it coincides with the (optimal) error term in ,fu.^^, . iii^i^^ii^i^, <j,o i^^j-i^ ^ ) 
the second term is just a higher order correction and the result is optimal as far as the order of the error 
term is concerned. However, for larger a the leading correction in (|3.9p is given by the second error term 
and it is due to the particular form of the trial function involved in the proof (see Section [4. 2p . 

In order to state a pointwise estimate analogous to p.7p , it is convenient to rescale the GP minimizer 
in the same way as when p^^ was obtained from p^^ by scaling. Thus we define (see also ()4.19p ') 



GP/ 



(3.10) 



with X = Ri^r. The scaled minimizer is concentrated in a neighborhood of a; = 1 and exponentially 
small everywhere else: 

Theorem 3.4 (Exponential smallness of the GP density, uj ^ 1/e) 
Set 



P = min 



4a(s + 2) 



, 1 



a{s + 2) 



3(s - 2) ' ' ■ s-2 
For any a, lui > , r £ = {f £ \ \1 — r\ > e^^^} and for e tending to zero, 



exp 



where 



1-^ 
3 



a{s + 2) 



(3.11) 



(3.12) 



(3.13) 



(Note that for both possible values of (3 the exponent 7 is positive.) Furthermore the density |\t'^^(a^)p 
converges in the sense of distributions to a Dirac delta function concentrated at x — 1. 

4 Proofs 

In this Section we prove the main results mentioned in Section [31 



4.1 The regime cu^e) ~ 1/e 

We start by proving the ground state energy asymptotics and the other results will follow as simple 
corollaries. The proof is quite similar to the proof of Theorem 2.1 in [CDYj : Like there, the main 
ingredient in the derivation of the upper bound for the energy is a trial function with a large number of 
vortices while the differences are essentially contained in a scaling argument. 

Proof of Theorem [37T] 

The lower bound is obtained again by simply neglecting the positive 'magnetic' kinetic energy in (|1.7p . 
namely 



> 



To get an upper bound we evaluate the GP functional on a trial function of the form 



(4.1) 



(4.2) 



where g^ is a phase factor, Xe a function that vanishes at the singularities of g^ and p^ a suitable 
regularization of p^^ . More precisely pe is defined as in Lemma 2.3 in |LSYlj . i.e., p^ = je * 



(4.3) 



Since ||je||i = 1, y/fh is L^— normaHzed. It is also clear that ps converges uniformly to p'^^ as e ^ and 
it is uniformly bounded in e, i.e., there exists a constant Cajg such that < C^^g. Furthermore, although 



Pe is not compactly supported, it is exponentially small in e for r sufficiently far from the support of p^^: 
For any f e R^, r > i?o 



Pe{r) 



•-out 1 

1 

2^ 



dr' exp 



r<-Rout 



1 



27re- 



■ exp 



r-Ro 



(4.4) 



Moreover, two different estimates for the gradient of hold true: By using the fact that |Vje| = e~^je, 
one can easily prove that \^ Pe\ < whereas by exploiting the regularity of p"^^ ■, i.e., ||Vyo"'"^|| < 

Ci^Q, one has ||V/3e||-^ < C^^g. Using both estimates we immediately get the bound 

||VV^||2 < C^Je. (4.5) 

The phase factor 5^ and the cutoff function Xe are defined as in [CDYj by placing vortices of degree 
1 at the points of the square lattice 



r < 2i?„„t - 2 



V24} 



(4.6) 



with spacing = Sy/e for some S > 0: Using complex notation ~ x + iy for r — (x, y) £ M'^ we define 



Cjec 

1 if \f-rj\>e'' 



Xe{r) = 



(4.7) 



(4.8) 



if \r — rj\ < e"^ 



for some 77 > 5/2. Note that the vortex lattice C has the same spacing as in |CDY) but it is extended to 
cover the whole of the support of p"^^ . The number, N^, of vortices and the normalization constant Cg 
satisfy the bounds < CujQ^s/e, due to (|4.6p . and 1 < < 1 + o{s^), since Xe < 1 and 77 > 5/2. By 
setting 

J6£ 



where Bi is a ball of radius e'' centered at r*, , we also have 



IVXell^ < 



1 



dr |V|f-f,||' < CN, < 



(4.9) 



The evaluation of the GP functional on gives 



2 £:TFr|^ 



< / dfp.Xe 



2 £TFn^ 



UJQ.S 



(4.10) 



where we have used the bounds (|4.5p and (|4.9p for the kinetic energies of and Xe- 

We can split the first term in (|4.10p into the contributions from B2B.aut and its complement respectively. 

Moreover, exploiting the pointwise estimate for r > 2i?out, 



< 



\^9e\ + ^<J2T^ + — 



< 



r3/2 



(4.11) 



and the exponential smallness (|4.4p , one has 



r>2Ro 



W -iAA g, 



< 



drr exp 



2fiou 



r-Ro 



< — ^exp 



Ro 



The remaining contribution can be estimated exactly as in [HDY) (see the proof of Theorem 2.1): 



\7-iAA g, 



lr<2Raut 

and an estimate similar to (|4.1ip yields, for f£ Bl 



+ C2 / drxi 



\/-iA,] g. 



< 



kec 



Ir-fll Ir-fjl infj^fe |fj - ffel Ir-fj] 



so that 



/ drxl (V - zA) 5/ < 2 / 



dfx?|V5e|' + 2 



dr 



A, 



^ Ci|Uj-s£g^'| C2N! ^ 



rl-2l) ' jr2-2)) 



The bound ()4.10p becomes then 

£^^[^e] <C I df 



^ df I (V - ile) 5e 



2 ^ £TF[|^^|2] c-.^,, 



(4.12) 



We now observe that the upper bound estimate for the first term in the r.h.s of the above expression can 
be simply taken over from Theorem 3.1 in [CDYj : A simple rescaling by 2i?out immediately yields 



dr 



- iA, ) g, 

and therefore, choosing 5 = 



/ df 


(v - iA, 


) 


/a 







^ < Cc:.o,A-|log£| 



(4.13) 



For the upper bound of the second term in (|4.12p we proceed as in [LSYlj : Denoting W{r) = —LOQr^ /A^ 
one has 

f TP [| VI/, 1 2] - £;TF < ^TF[^^] „ f TF[^TF] ^ < dr>TF(^) [^-^ ^ _ ^] (^) + ^(e) 

that is easily estimated using |iy(|f — £r'\) — W{r)\ < er'C{l + r^~^): 

[j,*W -W]{r) ^ ^ [ dr* [Wi\f-er^\)-Wir)]e-''' <eC'il + r''-^). (4.14) 
27r JR2 

We thus obtain the estimate f"""^ — E"^^ < CujqS and together with (|4.13p this finally yields the 
upper bound for the GP energy, i.e., s'^E^^ < e'^£'^'^[^e] < + C^^qEI loge|. 

□ 



Proof of Corollary \3A\ 

Defining 2a (r) = fi^^ — r'* + (ijjQ/4)r2 for all r > and using the negativity of a(r) outside the support 
of p^^ , we have 

/ df (iM/fPp - pTF)2 < f ^|^GP|4 _ 2a(r)|M/f + (p^^)^] = £TF[|vI/GP|2] _ ^tf 
since \\p^^\\l = f/^^ - E'^^ . The inequality ^^Fpcppj < e^^gGP g^j^^j Theorem O thus imply the result. 



□ 



Using Theorem 13.11 and Corollary 13.11 we can now show that the density of the minimizer is 
actually exponentially small outside the support of p^^ . 



Proof of Theorem [ 

The bound (|3.7p can be derived similarly to Proposition 2.4 in |CDY| or Proposition 2 in [AAB| . We 
present here only the proof of the first statement, since the second one is obtained exactly in the same 
way. 

The variational equation satisfied by '^f^ is 



where the GP chemical potential fi^ is fixed by jj^* 



GPlI _ 



1. Setting Us = and using 



Wo 

e 



*GP*^^GP 



. ,2 2|^GP|2 



4£2 



(4.15) 



we get 



The definition of /^e, Theorem 13.11 and CoroUarv 13.11 implv 



C^„e|loge|+ / dfl* 



,GP|4 



= + logel + / df - ip^^ f] < + a„ V^bi^. (4.16) 

On the other hand, since a'(i?out) < —Cloo < (see Eq. ()2.5p ). a simple Taylor expansion of a(r) in a 
neighborhood of _Rout yields 



for a possibly different constant > 0. Hence 



e^AU^ < 2 



(4.17) 



for any r € 6e = {f G | r > i?out +£^^'^/2} and e sufficiently small. Thus is subharmonic in 8^, so 
that, for any r and g with Bg{r) C 6^, 



Ueir) < / drUe < [ dfU^ 



1/2 



< 



drU'i 



r>Ra 



1/2 



< 



where we have used p.6p . If now we take r S 7^°"* and choose g — e^/"^/2, we have 



so that Ue{r) converges pointwise to in 7^°"*. Moreover from (|4.17p it follows that is a subsolution 
in r/"* of 

-Aw + C^^e~^l^w = 

(4.18) 

w;(5r-t) = C^„ei/6^1^^ 

whereas the r.h.s. of p.7p is a supersolution of the same problem for e sufficiently small. The result is 
then a consequence of the comparison principle. 



□ 



4.2 The regime Lj{e) > 1/e 



Proof of Theorem [ 

In order to capture the leading order term in the GP energy asymptotics it is convenient to rescale the GP 
functional in the following way: Setting ^'-'^[^I'] ee e^i?,^*^^^ [5'], with ^{x) = Rm^{r} and x = R^f, 
we have (remember that i?m depends on e through uja{e) — LUi/e") 



JM2 V / 



(4.19) 



where is the TF functional (|2.30p with ujq replaced with oji/e" and 

Ae = ^ 4 X X. 

The proof is thus similar to that of Proposition 12.11 By neglecting the (positive) first term in (|4.19p we 
get the lower bound 



(4.20) 



In order to obtain a corresponding upper bound we test the functional on a trial function "^^^e similar to 
the one used in the proof of Proposition l2.11 i.e., 



= yjkix) exp 



u;{s)Rl 



(4.21) 



where [ • ] stands for the integer part and the density p^{x) is defined in (I2.32p (we additionally require 
that ||Vv7||2 < oo). 

The estimate of the second term in (|4.19l) is thus already done in (|2.33p . It remains to bound the kinetic 
energy of i.e., 



u:ie)Rl 



[ dx (v - iAe) ^i,e = / dxVy^ + [ dx \- 
Smoothness of (and the assumption || V-\/j||2 < oo) yields the estimate 



ujie)Rl,x 



p^{x). 



dx 



while, for ^ sufficiently small. 



da?' 



Loie)Rl 



uj{e)Rl,x 



p^{x)<uj'{e)Rl / dx 



Jl-i/2 z 



dxh^< 



Altogether we get the bound 



(4.22) 



Optimizing with respect to the first two error terms we obtain the same error term as in (j2.31l) and the 
last term gives only a higher order correction, as long as a < ^^^^p-- On the other hand, for larger a, we 
consider the first and last terms in the above estimate and choose (in this case the second term can be 
neglected) 



which yields an overall remainder of order e e "-^ . 



□ 



Proof of Theorem [ 

The first part of Theorem 13.41 can be proved exactly as in the proof of Theorem 13.21 and we omit the 
details. The weak convergence to a Dirac delta function supported at a; = 1 is a simple consequence of 
the pointwise estimate p.l2p together with the L^— normalization of the density (see, e.g., the 

discussion in Section . 

□ 

5 Asymptotically homogeneous potentials 

We recall from the Introduction that a potential V{r) is called asymptotically homogeneous if there are 
constants k, c > such that the estimate 

\X-'V{Xr)-r'\<cX-''{l + r'') (5.1) 

holds for all A > 1 and all r E M+. We discuss here briefly how the results for the trapping potential 
can be extended to such potentials V{r) with suitable modifications of the error terms. 
The rescaling that produced (II. 2p leads in the general case to 

(5.2) 



i.e., the functional contains the rescaled external potential X^^V{Xr) with A = 

The estimates mentioned in Subsection 3.1 for cj <C 1/e generalize to asymptotically homogeneous 
potentials in exactly the same way as in [LSYlj . For the case uj ^ 1/e we have 

Proposition 5.1 (GP energy and density asymptotics for u ~ 1/e) 

Let the external potential V{r) > be asymptotically homogeneous of degree s > 2 in the sense of ()1.9p 
ana let i;GP and 'if^y denote the ground state energy and wave function of the functional (|5.2p . 
Then for any ujq, k > fixed and e tending to 0, 

= E'^^ + 0(e\ logel) + O (e^) . (5.3) 

Furthermore the density l^'^yP converges to p^^ strongly in L^{M.'^). 



Proof: The proof requires only a minor modification of the proof of Theorem 13. II Using ()1.9p . we can 
estimate 

e^\£^^[^]-£^^m<ce^ f df (1 + r^) |*|' , (5.4) 

so that the appropriate upper and lower bounds to E^y can be easily obtained: By testing the functional 
on 5"^^ we immediately get the upper bound E^y < E^^ + o(l), whereas taking "if = 5*^^ in the above 
inequality, one has the lower bound 



e^E'^^>e^Ef^-ce^- drr^\^%\\ (5.5) 
which yields the expected result, provided one can show that there exists a finite constant C^j^, such that 

dfrMvl/GPp <c^^^. (5.6) 



On the other hand, evaluating £y^ on a smooth radial function, we see that e^Efy < C', for some 



finite constant C' , so that 



but, using the trivial bound (ro = (cjq/2)=-2), 



r<ro 



drr' |*^^f < 



r>ro 



^0 



s-2 I 

+ 2 



together with p.9p . we get (|5.6p . i.e., 



^-0(1) dfr^l^rvl <CL„ + ^Vl^] +0(1). 



2 \ — 



The energy asymptotics follows then from Theorem 13. II 

In order to prove the ground state density convergence, it is sufficient to note that (|5.3p implies that 
jvP^yP is a minimizing sequence for the TF functional E"^^ . The statement can be thus obtained by a 

simple compactness argument together with identity of norms, = — 1 {see, e.g.. Theorem 

II.2 in |LSY2p . 



□ 



For uj{e) — uJi/e^^" we have the following generalization of Theorems 13.31 and 

Proposition 5.2 (GP energy and density asymptotics for lu ^ 1/e) 
Let the external potential V satisfy the same conditions as in Proposition \5.1l 
Then for any fixed uji, a, k > and e tending to zero, 



=-2 pGP _ 



i (l - + O (s^) + O (e s^) + O 



(5.7) 



Furthermore the rescaled density \'^^y{x)\ = R^\'^^y{Rmx)\ converges in the sense of distributions 
to a Dirac delta function concentrated at x = 1. 

Proof: It is sufficient to rescale the functional (jO)) as in (|4.19p . i.e., setting = e^i?,7/£y^[*], 

2 



da? 



R-'e—y [e^ — R^x^ - — 



and proceed as in the proof of Proposition 15. II to get the estimate 



e R„ 



Theorem 13.31 now yields the result. 

The rescaled density |\f^y(a?')| = i?^|5'^y(i?inS')| converges in the sense of distributions to a Dirac 
delta function concentrated at a: = 1 by the same arguments as before (cf. the proof of Theorem 13.41 in 
Section [ 



□ 



6 Conclusion 

We have analyzed in some detail the TF limit of the GP energy and density of a rapidly rotating Bose 
Einstein condensate in a two-dimensional trapping potential of the form r'*, s > 2. After discussing the 
scaling of the variables (that is necessary because of spreading due to the interaction and centrifugal 



forces) wc have estimated the energy with error terms whose order in the small parameter can be ex- 
pected to be optimal and proved the concentration of the GP density on the support of the TF density 
apart from exponentially small terms. The extension to asymptotically homogeneous potentials and the 
corresponding change of the error terms has also been discussed. 
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